REPORT  DOCUMENTATION  PA 

r990ft«A9  cO<i#a<OA  of  •nforfn«tion  >\  ntimatM  to  avt^a^a  l  ^our  OOf  '9 

gatharitW  matnia<A«n9  tha  data  waatd.  and  <0fnpfatin9  and  ravtawtng  tn#  <o4taaton  of  mf 
Sfction  of  infofnatton.  mciwdtng  »w9Mlton%  tof  r«ducin9  (hi$  Ourdan  to  ^aafttngton  Haadf 
Oavi«  Hid^nrav;  loita  t204,  Aritn9(on.  and  to  tna  OtfKt  o*  Mana^an^ant  and  tv 


AD-A256  645 

lilllilill 


I  data  vour<n 
a»oa<i  of  th.$ 
12 1)  >atfafM>n 


1.  AGIWCY  USI  ONLY  (tMw*  W#nit;  2.  RCPOMT  DATE  ...  _ _ _ 

July  1992  I  AIAA  Paper,  October  1991  to  July  1992 


4.  TITLE  AND  SUITITLE  5.  FUNDING  NUMBERS 

The  AlCl  Absorption  Feature  in  Solid  Rocket  Plume 
Radiation 


B.  AUTHOR($) 

McGregor,  W.  K. ,  Drakes,  J.  A.,  Beale,  K.  S.,  and 
Sherrell,  F,  G„ 


7.  PERFORMING  ORGANUATION  NAME($)  AND  ADOR£SS{E$) 

ArnoM  Engineering  Development  Center/DOT 
Air  Force  Materiel  Command 
Arnold  AFB,  TN  37389-5000 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AIAA-92-2917 


9.  SPONSORING /MONITORING  AGENCY  NAME(S)  AND  AOORESS(ES) 

U.  S.  Army  Research  Office 
P.  0.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


/^co 


11.  SUPPLEMENTARY  NOTES 

The  view,  opinions  and/or  findings  contained  In  this  report  are  those  of  the 
author (s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT  12b.  DISTRIBUTION  CODE 

Approved  for  public  release;  distribution  unlimited. 


13.  ABSTRACT  (Msximurri  200  words) 

A  phenomenon  observed  during  plume  diagnostic  measurements  acquired  during  simulated 
altitude  chcimber  tests  of  aluminum-loaded  solid-propellemt  rocket  motors  is  dis¬ 
cussed,  and  analytical  results  are  presented.  The  results  of  a  line-by-line 
calculation  of  the  aluminum  chloride  (AlCl)  feature  are  discussed  and  compared 
with  high-resolution  ultraviolet  CUV)  spectral  data. 


92-27619 

... ..... .Ml  iiiii  nil  Hi: 


14  SUBJEO  TERMS 

AlCl  absorption  Solid-Propellant  Rocket  Motor 

Plumes  Line-by-line  calculation 

Ultraviolet 


17.  SECURITY  CLASSIFICATION  IIS.  SECURITY  CLASSIFICATION  TlB.  SECURITY  CLASSIFICATION  I  20.  LIMITATION  OF  ABSTRACT 

OF  REPORT  I  OF  THIS  PAGE  I  OF  ABSTRACT  I 


IS.  NUMBER  OF  PAGES 


IB.  PRICE  CODE 


unclassified 


NSN  7540  01 -280-5500 


UNCLASSIFIED 


UNCLASSIFIED 


Standard  forin  298  (R*v  2-89) 
•>»k'i0M  a.  sia 

m  102 _ 


AIAA-92-2917 

The  Aid  Absorption  Feature  in 
Solid  Rocket  Plume  Radiation 

W.  K.  McGregor,  J.  A.  Drakes,  K.  S.  Beale,  and  F.  G.  Sherrell 
Sverdrup  Technology,  Inc.,  AEDC  Group 
Arnold  Air  Force  Base,  Tennessee 


■'iA  J.W 


'TED  1 


AIAA  27th  Thermophysics  Conference 

July  6-8,  1992  /  Nashville,  TN 


For  permission  to  copy  or  republish,  contact  the  American  Institute  of  Aeronautics  and  Astronautics 
370  L'Enfant  Promenade,  S.W.,  Washington,  D.C.  20024 
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Abstract 

An  absorption  feature  in  the  mid-ultraviolet 
spectral  region  of  the  radiation  from  aluminized 
solid-propellent  rocket  plumes  has  been  identi¬ 
fied  as  the  X'E  -*  A’n  transition  of  aluminum 
chloride  (AlCI).  A  proposed  model  of  the  radia¬ 
tive  transfer  has  the  continuum  radiation  from 
the  aluminum  oxide  particulate  core  passing 
through  the  outer  gas  boundary.  A  line-by-line 
calculation  is  then  performed  in  a  temperature- 
density-path  length  parameter  space  to  attempt 
a  tit  to  the  data.  The  final  temperature,  density, 
and  path  length  agreed  well  with  flow-field  calcu¬ 
lations.  The  spectral  agreement  between  model 
and  data  was  marginal,  but  sufficient  spectral 
structure  and  band  shape  fits  were  found  to 
merit  further  work  on  this  mechanism  as  a 
diagnostic  tool  in  solid-propellant  rocket  studies. 


Introduction 

An  interesting  feature  in  the  mid-ultraviolet 
spectrum  of  the  radiation  from  aluminized  solid- 
propellent  rocket  exhausts  has  been  noted 
several  times  in  unpublished  reports.  This  feature 
is  a  deep  absorption  near  262  nm  of  the  other¬ 
wise  continuum  spectrum  typical  of  hot  aluminum 
oxide  particles.  It  has  been  attributed  to  the  aluminum 
chloride  (AlCI)  A-X  transition,  although  the  spectral 
resolution  was  not  sufficient  to  resolve  bandheads.’ 
Recently,  spectra  were  obtained  having  sufficient 
resolution  to  identify  the  vibrational  bands  and  these 
will  be  reported  here.  Curiously,  spectra  obtained 
from  an  on-board  spectrometer  viewing  a  rocket 
plume  having  an  aluminized  propellent  also  contained 
this  feature. 2 

One  obvious  question  arising  from  such  an 
inherent  spectral  feature  is  “Can  it  be  used  to  diag¬ 
nose  any  properties  of  the  rocket  exhaust  flow?”  To 


that  end,  a  modeling  effort  was  undertaken  and  this 
paper  reports  on  the  success  of  that  work. 


Fig.  1 .  Schematic  of  solid  rocket  motor  test  configuration 
used  for  UV  spectral  measurements. 


Description  of  the  Spectra 

When  a  moderate  resolution  (—0.8  nm)  ultra¬ 
violet  scanning  monochromator  is  pointed  normal  to 
the  axis  of  an  exhaust  plume  from  an  aluminized 
solid-propellant  rocket  motor  being  fired  in  an  altitude 
test  cell  (Fig.  1),  the  spectra  shown  in  Fig.  2  results. 
The  continuum  on  either  side  of  the  feature  at  262 
nm  and  the  thermal  behavior  of  the  radiation  in  the 
visible  and  infrared  regions  suggests  that  the  feature 
is  due  to  absorption  of  a  gaseous  specie.  A  quick 
examination  of  the  predicted  products  of  combustion 
of  ammonium  perchlorate  with  aluminum  and  hydro¬ 
carbon  solids  along  with  a  trip  through  Pearce  and 
Gaydon’s3  atlas  of  spectra  reveals  that  the  most 
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Fig.  2.  Moderate  resolution  mid-UV  spectra  from 
solid-propellant  rocket  motor  plume. 

likely  absorber  is  aluminum  chloride  (AlCI).  To  investi¬ 
gate  this  identification  more  thoroughly,  higher-reso¬ 
lution  spectra  were  obtained  during  a  subsequent 
test.  Because  of  the  small  vibrational  (~450  cm-’) 
and  rotational  (~0.25  cm-’)  constants  of  both  AlCI 
states,  a  spectral  resolution  of  about  0.1  nm  was 
required  to  produce  the  band  structure.  The  resulting 
spectrum  (Fig.  3),  indeed  does  show  the  expected 
structure,  and  many  of  the  bands  listed  in  Pearce 
and  Gaydon,  Table  1,  are  identified. 


WAVELENGTH,  nm 

Fig.  3.  High  resolution  UV  spectra  from  solid-pro¬ 
pellant  rocket  motor  plume. 


The  AlCI  Molecule 

Aluminum  chloride  sublimates  at  temperatures 
above  460  K.  A  calculation  of  the  plume  properties 
using  the  JANNAF  Standard  Plume  flow-field  codes^ 
gives  the  solid  and  gaseous  temperature  profiles  of 
Fig.  4  for  this  particular  rocket  plume.  The  minimum 
predicted  temperature  is  about  600  K  so  that  the  AlCI 
IS  gaseous.  Furthermore,  the  predicted  mole  fraction 
of  AlCI  IS  about  2.6  percent,  and  should  follow  the 
static  pressure  radial  profile  (Fig.  4). 


Table  1.  Bandheads  of  the  AlCI  A’n  -  X’E  System 
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NOTE:  THE  LETTERS  R  AND  V  AFTER  THE  WAVELENGTH  INDIUTE 
THE  DIREQION  OF  DEGRADATION  OF  THE  BAND,  WHEREAS 
THE  NATURE  OF  THE  HEAD  (R  AND  Q)  IS  INDIUTED  BEFORE 
THE  VIBRATIONAL  QUANTUM  NUMBERS. 

REF  :  JHimmmaimofmiKuiAitsPinu.i.ii  b  pearse 
AND  A.  G.  GAVOEN,  CHAPMAN  AND  HALL  LTD.,  1950,  p.  47. 


Fig.  4.  Predicted  radial  profiles  of  temperature  and 
pressure  for  solid-propellant  rocket  motor  at 
axial  location  where  UV  spectra  were 
measured  (standard  plume  flow-field  code 
using  three  particle  groups). 
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The  molecular  constants  of  the  two  electronic 
states  (A'n  and  X’S),  measured  by  several  inves¬ 
tigators  and  reported  critically  by  Huber  and 
HerzbergS,  are  given  in  Table  2.  The  values  are  for 
AI35CI,  although  AI37CI  also  naturally  occurs;  these 
two  isotopes  occur  in  the  ratio  75.5:24.5.  The 
relatively  small  vibrational  constant  of  the  ground 
state  (480  cm-’)  leads  to  a  vibrational-rotational 
fundamental  at  a  wavelength  of  about  21  pm. 
(Observation  of  this  LWIR  emission  band  from  rocket 
plumes  has  not  been  reported.)  The  constants  of  the 
X’S  state  have  been  re-measured  recently  in  the 
infrared  bands,  and  improvements  m  accuracy  have 
been  reported.  These  data  will  be  introduced  subse¬ 
quently.  The  rotational  constants  for  this  closely 
packed  molecule  are  also  quite  small  (-0.254  cm-’), 
and  the  overlapping  branches  make  resolution  of  the 
lines  near  impossible. 

Table  2.  Molecular  Constants  for  the  Ground  and 
Resonant  States  of  APSQ  (cm-’) 


PLUME  CROSS  SEQION 


Fig.  5.  Physical  model  of  solid-propellant  rocket 
motor  plume  UV  spectral  measurement 
experiment. 

A  complete  solution  to  the  monochromatic 
radiative  transfer  equation  is  desired.  The  equation  is 
written 
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The  vibrational  correction  to  the  rotational  constant 
given  as  footnotes  in  Ref.  5  is  a  necessity  for  proper 
modeling  of  this  band.  Also,  because  of  the  small 
values  of  Bg,  rotational  quantum  numbers  out  to  at 
least  100  must  be  used  at  the  temperature  of 
interest. 

Radiative  Transfer  Model 

The  preliminary  model  adopted  for  this  study  is 
illustrated  m  Fig.  5.  It  is  supposed  that  the  core  of 
solid  particulates  provides  a  continuum  source  for 
transmission  through  a  layer  of  gas  containing  the 
absorbing  AlCI.  The  continuum  radiation  from  the 
particulates  was  curve  fit  by  joining  the  spectra  at 
both  extremes  of  the  absorption  feature,  and  served 
as  a  source  for  the  absorption  calculations.  The  gas 
layer  was  considered  uniform  for  these  preliminary 
calculations,  with  temperature  and  number  density  as 
parameters.  The  absorption  spectra  were  modeled  as 
Voigt  lines  with  the  spectrometer  triangular  slit  func¬ 
tion  considered  to  be  much  greater  than  a  line  width. 


where 

Ntot  is  the  total  number  density  of  emitter/ 
absorber. 

Lv  is  the  line  shape  function 
Se  is  the  molecular  emission  strength 
Sa  is  the  molecular  absorption  strength 
1^  is  the  spectral  radiance,  (i.e.,  W/sr-cm2- 
cm-’)  at  wavenumber  u 

X  is  the  spatial  coordinate  within  the  media 

The  boundary  condition  at  x  =  0  is  given  in  this  case 
by  the  fit  to  the  continuum  spectra,  neglecting  the 
AlCI  absorption  feature;  i.e.,  the  continuum  back¬ 
ground. 

Conceptually,  the  calculation  involves  3  steps: 

1.  Determine  the  positions  of  all  spectral  lines. 
Determine  the  AlCI  emission  anu  absorption 
strengths  for  each  line,  based  on  a  rotational  and 
vibrational  temperature.  This  yield';  a  spectral  line 
atlas  of  line  strengths  versus  line  center  wavenumber 
for  all  lines. 

2.  Given  Sa  and  Se-  and  Nfor.  solve  Eq.  (1). 
Here,  only  one  spatial  zone  is  considered  so  that 
integration  of  Eq.  (1)  is  simplified. 

3.  Convolve  the  result  of  2  for  each  line  with 
spectrometer  response  function.  In  this  case  a 
triangular  slit  function  suffices  since  the  actual  line 
width  IS  much  smaller  than  the  instrument  mono¬ 
chromatic  response  function. 
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Part  1 

The  strength  of  each  emission  line  is  given  by: 


-  A  TT' 

E~  4n  vv  2J'  +  l 


(2) 


x{A;v‘;J') 


where 

A 


v  V  are  the  Einstein  coefficients^ 

Sj'j"  are  Honl-Londons8 

V  wavenumber 

x{A;v’;J’)  is  the  fractional  population  of  the 
A’n  state  for  v‘  and  J'  level,  nor¬ 
malized  by  the  partition  function,  i.e. 
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Table  3. 

Einstein  Coefficients  for  the 

Ain 

-  X^  + 

transition  in  AlCl  (Ref.  7) 
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VALUES  NOT  GIVEN  ARE  LESS  THAN  0.5  PERCENT  OF  VALUE  AT  (0,0) 


energy  of  quantum.  Now,  in  terms  of  the  vibrational 
quantum  number  v,  the  rotational  quantum  number  J, 
and  the  electronic  angular  momentum  quantum 
number  A, 

E,,  =(ji  (v  -I-  —  a>  X  (v  -f-  -t-  w  y  (v 

V.  6  6*^6 


where  g,  is  the  degeneracy  factor  of  the  given  state. 

E„  =B  [JfJ-t-D- A^l- D  [(J(J  +  1)-AT] 

rt.  V  V 

The  absorption  strength  is  similarly  given  by  ' 
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with  the  Einstein  B  coefficient  as 
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(5) 


The  band  strength  (or  lifetime)  of  the  AlCI  A’n  -► 
X'S  (0,0)  transition  has  been  measured®  to  be  6.4  x 
10^  sec, where  no  effort  was  made  to  separate  the 
two  isotopic  states.  The  relative  band  strengths  were 
calculated  by  Langhoff  and  Bauschlicher^  from  first 
principles.  They  are  given  as  Einstein  coefficients  in 
Table  3.  Note  the  relatively  strong  values  for  the 
Av  =  0  sequence,  indicating  that  these  will  be  the 
most  evident  in  the  structure,  as  also  reported  in  the 
experimental  analyses. 3  These  calculations  of  band 
transition  probabilities  were  made  for  AI35CI  with  the 
assumption  that  the  results  are  the  same  for  AI37CI. 

Energies  of  individual  states  are  given  by 

E  =  T  -1-  E„  -f-  E„  (6) 

11  V  R. 


D  =  D  B  (v  -I-  A) 

V  8  ~  8  ^ 

For  mass  A)37Cl  all  constants  are  modified  by  the 
ratio  p  =  V35/37 


g(37)_g(35)  2 
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(37)  (37) 

a  =a 
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3  . 
p  etc. 


The  molecular  constants  for  the  A’ll  state  were 
taken  from  Huber  and  HertzbergS  and  are  listed  in 
Table  2.  For  the  ground  electronic  state  the  latest 
reported  measurementss  use  the  Dunham  expres¬ 
sions  for  the  sum  of  the  vibrational  and  rotational 
energy. 


T^.j=  X 

l,m 


(9) 


where  T,  is  electronic  state  energy,  Ev  is  the  ^he  values  of  Y,^  were  given  to  the  author  by 

vibrational  state  energy,  and  Er,  ,s  the  Rotational  Bernaths  and  appear  m  Table  4. 
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Table  4.  Constants  for  Dunham  Level  Expression 
of  Energy  Levels  of  the  state  of  AlCI 


Moleculai  Constants  for  AICI(X'!:)  used  in  Term  Energy  Expression: 
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VALUES  ARE  IN  UNITS  OF  CM-'  (REF  8) 


With  the  transition  energies  defined,  line  positions  are 
obtained  by  subtraction 

=  E{A,v',J')  -  E(X,  v",J") 

A  line  atlas  is  built  up  from  these  expression  for  v', 
v"  out  to  8  and  J‘,  J"  out  to  200  for  both  isotopes. 

Part  2 

The  line-by-line  model  then  solves  Eq.  (1)  tor  each 
line  in  the  atlas.  A  Voigt  profile  for  each  line  is 
assumed.  The  Se  and  Sa  are  calculated  at  a  mini¬ 
mum  of  ten  spectral  frequencies  between  spectral 
line  centers  by  summing  contributions  of  all  relevant 
lines  at  that  position.  Here  v’.v”  out  to  8  and  J'  J" 
out  to  200  are  used  to  generate  a  large  number  of 
stored  numbers.  Thus,  for  20.000  spectral  lines,  at 
least  200,000  radiances  and  transmittances  are 
generated.  The  program  also  does  spatial  integration, 
although  for  this  use  (1  zone),  this  is  trivial. 

Part  3 

The  spectrometer  is  modeled  by  a  triangle  function  of 
2A  full  width  at  half  maximum.  The  convolution  of  the 
integral 

Ki-- i'')G(Odi'’ 


is  done  numerically  using  the  trapezoid  rule,  where 
Hfv)  represents  the  instrument  reading  at  »'  and  G(;'‘) 
represents  the  instrument  function, 

Results  of  Modeling 

The  approach  taken  was  to  vary  the  vibrational 
and  rotational  temperatures,  the  column  density,  and 
the  column  length  in  Eq.  (2)  until  the  closest  fit  to  the 
experimental  data  was  achieved.  The  temperature 
and  pressure  were  varied  about  the  predicted 
averages  along  the  column,  (see  Fig.  4)  and  the  pre¬ 
dicted  mole  fraction  of  AlCI  was  used  initially.  It  was 
first  necessary  to  calculate  the  constants  accurately 
for  each  vibrational  state  of  each  of  the  isotopic 
molecules.  This  proved  to  be  very  sensitive  in  obtain¬ 
ing  a  fit  to  the  vibrational  bandheads.  Then,  the 
number  of  J  values  necessary  to  fit  the  absorption 
level  was  determined. 
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Fig.  6.  Comparison  of  computed  with  measured 
spectra  of  SRM  plume  radiation  at  predicted 
plume  conditions. 

Only  two  results  will  be  shown:  First,  the  result 
using  the  predicted  temperature  and  AlCI  density. 
Fig.  6.  and  the  final  result,  the  best  fit.  Fig.  7. 
Referring  to  Figs.  4  and  5,  the  parameters  used  for 
the  baseline  calculations  were 

Vibrational  Temperature,  Ty  =  800  K 
Rotational  Temperature,  T,  =  800  K 
Number  Density,  ri  =  2.36  x  10’6  cm-3 
Absorber  Path  Length,  ^  =  5  cm 

Quite  clearly  from  Fig.  6,  the  fit  of  the  calculation 
to  the  data  is  poor.  The  structure  is  much  smoother 
than  the  data,  although  the  same  spectral  resolution 
was  used  m  the  model  as  for  the  spectrometer.  In 
the  process  of  iterating  on  temperature  and  density. 
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Fig.  7.  Comparison  of  computed  with  measured 
spectra  of  SRM  plume  radiation  at  "best  fit" 
conditions. 

a  better  fit  could  be  found  m  the  central  portion  of  the 
absorption,  but  the  structure  in  the  wings  became 
poorer.  It  was  this  lack  of  spectral  fit  that  led  to  use 
of  the  more  accurate  data  for  the  X'S  state.  Finally, 
since  the  molecular  data  for  the  upper  state  are 
somewhat  dated  and  some  uncertainty  exists  in  the 
relative  values  of  the  rotational  constant  for  the  two 
isotopes,  the  vibrational  and  rotational  constants  of 
the  A‘ri  state  were  also  iterated.  Since  this  process 
could  go  on  and  on.  a  final  fit,  Fig.  7,  was  produced. 
The  final  parameters  are  as  follows: 

Vibrational  Temperature,  Ty  =  1,000  K 
Rotational  Temperature,  Tr  =  1,000  K 
Number  Density,  i]  =  2.36  x  lO’S  cm-3 
Absorber  Path  Length,  ^  =  5  cm 

and  the  molecular  constants  were  also  shifted 

u)e  450  -*  430  cm  1 
WeXe  4.37  -»4.1  cm  ' 

Be  0.259  -»0.249  cm  ’ 

Here,  the  spectral  match  is  fair  and  the  overall 
absorption  matches  pretty  well,  so  that  continued 
Iteration  appears  unwarranted. 

The  final  result  is  reasonable.  Because  of  the 
high  particle  temperature,  it  is  not  surprising  to  have 
a  gas  temperature  somewhat  larger  than  predicted 
(1,000  K  versus  800  K  predicted).  Also,  because  of 
small  particles  in  the  outer  gas  flow  it  is  not  sur¬ 
prising  that  the  radiation  from  the  gas  layer  masks 
the  absorption  at  the  deepest  absorption  frequency. 
One  might  expect  that  an  arbitrary  shift  of  the 
wavelength  for  the  experimental  data  would  produce 
a  better  fit.  However,  attempts  to  make  this  match 
were  only  marginally  successful.  The  simple  physical 


model  of  Fig.  5  cannot  be  correct,  since  the  radiation 
is  obviously  extended  and  the  AlCI  exists  all  through 
the  plume.  However,  the  modeling  brought  assurance 
that  the  absorber  was  AlCI,  and  only  a  slight  variation 
of  temperature  brought  the  model  results  very  near 
the  measured  data. 

Concluding  Remark 

Clearly,  if  the  AlCI  feature  is  to  be  used  for 
diagnostics,  laboratory  work  using  a  heated  absorp¬ 
tion  cell  will  be  required  to  reach  a  better  fit  between 
calculated  and  measured  spectra.  In  a  heated  cell 
experiment  the  temperature  and  pressure  parameters 
would  be  well  known,  leaving  only  the  spectrometer 
wavelength  calibration  and  the  molecular  constants 
unknown.  The  issues  raised  above  can  then  be 
resolved. 
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